Non-equilibrium Scaling Properties of a Double Quantum Dot System: Comparison 
between Perturbative Renormalization Group and Flow Equation Approach 
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Since the experimental realization of Kondo physics in quantum dots, its far-from-equilibrium 
properties have generated considerable theoretical interest. This is due to the interesting interplay 
of non-equilibrium physics and correlation effects in this model, which has by now been analyzed us- 
ing several new theoretical methods that generalize renormalization techniques to non-equilibrium 
situations. While very good agreement between these methods has been found for the spin- 1/2 
Kondo model, it is desirable to have a better understanding of their applicability for more compli- 
cated impurity models. In this paper the differences and commons between two such approaches, 
namely the flow equation method out of equilibrium and the frequency-dependent poor man's scal- 
ing approach are presented for the non-equilibrium double quantum dot system. This will turn out 
to be a particularly suitable testing ground while being experimentally interesting in its own right. 
An outlook is given on the quantum critical behavior of the double quantum dot system and its 
accessibility with the two methods. 

PACS numbers: 05.10.Cc, 73.63.Kv, 72.10.Fk, 75.30.Hx, 72.15.Qm 
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I. INTRODUCTION 

The recent advances in nanotechnology permit to 
probe far-from-equilibrium transport properties of cor- 
related electron systems. The paradigm for such exper- 
iments is the realization of Kondo physics in Coulomb 
blockade quantum dotsi. This combination of non- 
equilibrium with correlation effects is theoretically chal- 
lenging and has led to considerable efforts to develop suit- 
able theoretical tools. The key challenge is that beyond 
linear response theory one does in general not know how 
to systematically construct the steady current-carrying 
state, except by following a difficult real time evolution 
problem. 

During the past five years new powerful methods, both 
numerical and analytical, have been developed and used 
for investigating non-equilibrium quantum many-body 
problems. As analytical renormalization group methods 
have played an enormous role in understanding corre- 
lation physics in equilibrium, a number of these meth- 
ods are directly based on renormalization ideas and carry 
them over to non-equilibrium situations: the frequency- 
dependent renormalization group 2 ' 3 , the real time renor- 
malization group^ and the flow equation method^—. 
While these methods show very good agreement for the 
Kondo model, it is clearly desirable for future applica- 
tions to understand their relation and respective advan- 
tages in more detail. 

In this paper we address this question for the 
frequency-dependent renormalization group and the flow 
equation method in the case of a more complicated quan- 



tum impurity model, namely for a double-dot system in 
the Kondo regime. In this system two Kondo dots are 
coupled via a spin-spin interaction, which is a setup re- 
lated to recent experiments^^. In addition, this model 
is particularly interesting in the present context for two 
additional reasons: i) It is known to exhibit an inter- 
esting quantum phase transition in equilibrium— and 
one expects non-equilibrium properties to be especially 
important at quantum phase transitions, ii) Based on 
the investigation of the non-equilibrium spin- 1/2 Kondo 
model, one knows that the decoherence generated by the 
steady state current plays the key role in understanding 
the far-from-equilibrium properties 2 -^. Now decoher- 
ence enters via two seemingly very different mechanisms 
in the frequency-dependent renormalization group and 
the flow equation method. For the frequency-dependent 
renormalization group one identifies suitable Korringa- 
like decoherence rates that are then used to explicitly cut 
off the renormalization flow^i^. On the other hand, in 
the flow equation approach decoherence terms appear as 
two-loop contributions in the scaling equations^ 7 -. For 
the conventional spin- 1/2 Kondo system the different 
renormalization-based approaches show very good quan- 
titative agreement including line-shapes and Korringa 
rates^i^. 

The difference how decoherence enters highlights 
the respective advantages of these approaches: the 
frequency-dependent renormalization group (like real 
time RG) has a straightforward diagrammatic represen- 
tation, while the flow equation method treats both many- 
particle coherence effects (like Kondo physics) and deco- 
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herence on the same footing in the scaling equations. In 
the following we will analyze how these differences man- 
ifest themselves in the non-equilibrium double-dot sys- 
tem. Previous studies of the non-equilibrium double-dot 
system based on the frequency-dependent renormaliza- 
tion group and non-equilibrium perturbation theory have 
been published in Refs. Il^filil . 

This paper is organized as follows. In section UU we 
first introduce the Hamiltonian of the double quantum 
dot system and discuss the various approximations used 
in its derivation. Section Mil explains the flow equation 
method and shows its application to the double dot sys- 
tem. In section IIVI we then introduce the perturbativc 
renormalization group approach and discuss the scaling 
equations for our system. After comparing the two meth- 
ods in leading logarithmic order in section [Vl we discuss 
generalizations of the two approaches in section IVTI which 
include decoherence effects. At the end of section IVII 
we then compare the two approaches again and show 
that the results at the decoherence scale are in very good 
agreement, although the underlying methodologies are 
quite different. 



II. MODEL 

The simplest model to illustrate the competition be- 
tween a spin singlet and a Kondo singlet formation is 
the two- impurity Kondo model. It has been studied in 
detail Ref. [Xl| and gained new life with the progress in 
nanotechnology and the possibility to study two single 
quantum dots interacting with each other— i 10 i 16 . 

Here we study the problem of two quantum dots where 
two artificial impurities are attached to leads, in contrast 
to the historical two-impurity model which contains two 
impurities embedded in a metalii. 

The double quantum dot (DQD) model describes two 
spin- 1/2 degrees of freedom denoted as Sl and Sr, 
which are each Kondo coupled to conduction band elec- 
trons with an additional mutual spin exchange interac- 
tion H CK — K SlSr, which is assumed to be antiferro- 
magnetic K > 0. The Hamiltonian of the system is given 
by 

H = H i oa d s + H cx + -ff Kondo (1) 

The conduction band electrons are described by 

Pleads Yl £k >J : C L>%<7 » ( 2 ) 

j k,a 

where the lead index j is specified later on, e^j is the en- 
ergy of the corresponding conduction band electron and 
c kjc ( c kja) are the corresponding creation (annihilation) 
operators for a conduction electron with momentum k 
and spin a. The notation : . . . : denotes normal ordering 
with respect to the non-interacting ground state. 



The Kondo interaction with the leads is 

#Kondo J k'i : S LS(k'j)(kj) ■ 

j k'k 

+ 22z2 J k>k '■ Su.S{k'j)(kj) ■ ■ (3) 

j k'k 

where J^,\ (J^ k ) is the coupling of the left (right) quan- 
tum dot spin to the spin density of the conduction band 
electrons in the lead j 

■ S(k'j)(kj) ■ = \foo' ■ c\'ja' C kja 1 ( 4 ) 

a' a 

and r are the Pauli matrices of a spin-1/2. 

In the two- impurity model the spin-spin interaction be- 
tween the Kondo spins is mediated by the RKKY inter- 
action^, i.e. for antiferromagnetic coupling generated by 
two Kondo spin-flip interactions. The RKKY interaction 
depends in sign and strength on the distance between 
the two impurities and since it is an effective interac- 
tion in J 2 it is both retarded and small (at least every 
reasonable theory should do so). Note that the RKKY 
interaction in this case scales to the same degree as the 
Kondo interaction in scaling theory. 

On the contrary singlet-triplet states in quantum dots 
can arise from other physical effects, for example from 
exchange couplings and/or orbital degeneracies. There- 
fore the effective spin-spin interaction between the Kondo 
impurities can be tuned independently from the Kondo 
interaction with the leads^. 

For both methods we therefore include the spin ex- 
change interaction in the unperturbed Hamiltonian 

Ho =#lcads + -ffcx, (5) 

and we treat the Kondo interaction -f/Kondo as a small 
perturbation. The eigenstates of the unperturbed double 
dot spin system are singlet | 0, 0) with a total spin S = 
and triplet states | 1, to) with a total spin of S — 1 where 
to = { — 1,0,1}. The perturbative RG focuses on the 
flow of a generalized coupling function and the scaling 
equation is derived by diagrammatic perturbation the- 
ory in the vertex. In order to do the perturbation theory 
a pseudo-particle representation for the spin operators 
is introduced. In the flow equation method the scaling 
equations are derived from infinitesimal unitary trans- 
formations. These involve mainly the commutation rela- 
tions of operators and therefore the spin can be treated 
as an operator. 



A. Pseudoparticle representation 

In order to calculate diagrams in perturbation theory 
and also in perturbative RG, we introduce the pseudo 
particles dl which create a singlet or triplet state 7 € 
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{s, f_, to, t + }. The spin exchange interaction Hamilto- 
nian is thus diagonal in the pseudo particle operators, 



-Kd\d s 
4 



-K V<4 d t 

7 



(0) 



The left and right spin (upper and lower sign, respec- 
tively) can be represented by bond operators 1 ^: 



S L/R - 2^ d t+ dt + 



9+ 



(S 



L/RJ 



1(4 



t d t . 



- d\_d t _ ±d\d to ±d\ d s ), 
(7) 

■d] + d t0 ±dtd t _Tdt + d s ). 

(8) 



The constraint 
Q = d\d s 



d\ + d t+ 



4A 



d\_d t _ 



1. 



(9) 



is fulfilled by calculating the physical observable from the 
expectation value 



(0)q= 



lim 



(QOh 
(Q)x 



(10) 



where (■ ■ -)\ is the average over the Hamiltonian H\ = 
H + XQ where the constraint enters as a chemical poten- 
tial which is set to infinity at the end of the calculation™,. 

It is convenient to introduce a matrix representation 
of the Kondo spins in the bond operator notation by 
defining a generalized Pauli matrix T a leading to 



= E ^4' (^) 7 ' 7 d 7 



(11) 



In the case of the exchange coupled double quantum dot 
system this generalized Pauli matrices are given by 



/0 1 
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(12) 
(13) 

(14) 



Note that the lower right 3x3 matrix, i.e. the triplet 
states, represents the Pauli matrices for a spin-1. 

Using this notation the interaction Hamiltonian Q is 
given in the general form 



-Hint = 



^ V 4 

a,j—L : R 7 7 k' a' ;ka 



7'7 



x di,d^ : r 



7 7 



k' ja' kja 



(15) 



During the renormalization we will find that the cou- 
pling between triplet states J u d\ m d t , flows differently 
than for the spin couplings including a singlet-to-triplct 
transition J s td\ d s and J ts d\d t . 



B. Spin notation 

If we do not introduce pseudoparticle states, but keep 
the spin operator as a quantity, we find that transitions 
between the eigenstates of the DQD are given by 



(s L + Sr) \S,m) -> IS 1 , m), 



S L + S R )\S,m)->\S,m), (16) 
(S'l - Sr) + 2i(S L x S R )j |0, 0) -> |1, to), (17) 
(S L - S R ) - 2i{S L x Sr)) |l,m) -> |0,0). (18) 



Therefore we rewrite the Hamiltonian in Eq. ([3]) for 
the flow equation treatment by 

Hurt = E E '■ + g {k']){kj) ■ 

j=L,R k'k 

+ E E P L ■■ ((Sl Sr) + 2 i (S L x S R )) s [k , j)m : 

j=L,R k'k 

+ E E M l'k ■■ ((Sl Sr) - 2i(S L x S R )) s (k , Mkj} :, 



j=L,R k'k 

where the couplings are defined by 

Tsum,j _ \ ( T L,j . T R,j\ 
J k'k — 2 \ k'k ^ J k'k J ' 



pi — 1 | I ( i 1 • _ i" • i ... c i 

r k'k — 2 \ 2 v ' fc ' 



J 

k'k I ' 



M « = \ (\ 



(4 



"k'k 



1 ■)!',') Qi, 



k'k I ' 



(19) 

(20) 
(21) 
(22) 



The interaction Q k , k : 2i{Sl x Sii)i?(k'j)(k3) '■ i g P er se 
not present in the initial setup, but it turns out that the 
system of equations does not close if it is not included. 
This leads to the following initial conditions for a general 
flow parameter B 



J L k'i(B 


= 0) 


— J k'k 


(23) 


4i(B 


= 0) 


- J R ' j 
~ J k'k 


(24) 


Qi'kiB 


= 0) 


= 


(25) 



The following symmetry relations have to be fulfilled 
during the flow due to the hermiticity of the Hamiltonian 



J k'k ~ J kk' ' 

pL = Kk- 



(26) 
(27) 



The interaction P k , k refers to a scattering process in- 
volving a singlet to triplet transition as indicated in 
Eq. dT71) . The interaction M 3 k , k is the hermitian con- 
jugate of Pl k = (M{ lk )\ 
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C. Discussion of the leads 

In this paper we will concentrate on the case of two 
exchange coupled quantum dots which are not coupled 
symmetrically to a set of leads. We want to test the two 
quantum dots independently and therefore we assume 
that there are two leads attached to each quantum dot 
such that transport can take place through each quantum 
dot independently, see Fig. [TJ 




FIG. 1: Double Quantum Dot Setup: the residing electrons 
with a spin degree of freedom on the left (L) and right (R) 
quantum dot are coupled mutually by the exchange interac- 
tion K. Two lead 1 and 2 (3 and 4) measure the current 
through the left (right) quantum dot. The leads are coupled 
to the quantum dots only by the Kondo spin exchange inter- 
action. In the flow equation method we symmetrize the leads 
1 and 2 (3 and 4) to one symmetric left (right) lead as denoted 
by dashed lines. 

A completely symmetric setup corresponds to two 
Kondo impurities embedded in a metal, but in this case 
the two impurities couple to the same position in the lead 
and for example the RKKY interaction is not defined 
properly. On the other hand for completely symmetric 
coupling the initial conditions are J^™ J (-B = 0) = Jo 
and P 3 k , k {B = 0) = M J k , k (B = 0) = 0. The singlet-triplet 
transitions are not allowed in the beginning and never 
created during the flow. Note that [(S L + S R ) 2 ,H] = 
and thus the total spin is conserved. The flow yields the 
physics of a non-degenerate singlet or a spin-1 coupled to 
leads depending on the initial condition. 

In order to make non-equilibrium electron transport 
possible, two leads have to be put at different chemical 
potential, fii/2 = ±eVL/2 where I and 2 denote the two 
leads attached to the left dot; similarly 3 and 4 for the 
right two leads, ^3/4 = ±cVr/2. While we treat the 
four leads independently in the pRG approach, we can 
introduce a simplification due to symmetry arguments in 
the flow equation calculation. 

We focus now for a short while on the left dot and 
leads 1 and 2. It has been shownii, that a quantum dot 
is coupled only to the even mode of the two leads defined 
as 

cl = Jt^-ci + Jri— 02 ( 28 ) 
V I + r L V l + 

where — J11/J22 (tr = J33/J44) is the asymmetry 
parameter. For simplicity we only discuss rt = r R = 1 



in the numerics. Note that the extension to r a 7^ 1 is 
trivial. Therefore we have to take into account only 2 
instead of 4 leads but with a double step-like occupation 
function n a (e) 

( if e > eV a /2 

n a {e) = I r a /{l + r a ) if |e| < eV a /2 . (29) 
[ I if e < -eV a /2 

where a — L,R denotes the left or right set of leads. 
A non-zero voltage enters the calculation thus via the 
normal ordering of the lead electrons, see Ref. @. In the 
following we assume that initially no cross-talk between 
the left dot and the right leads and vice versa is present, 
such that we can "drive" and "probe" the two quantum 
dots independently by each their leads. This corresponds 
to the initial conditions 

■#fc(B = 0) = Jjg£(B = 0) = 0. (30) 

Consequently the initial conditions read 

jft(B = 0) = jjff(B = 0) = 2J , (31) 
=* 4T L (B = 0) = 4™> H (B = 0) = J Q , (32) 

P£, k (B = 0) = M% k (B = 0) = ± Jo , (33) 

P$ k (B = 0) = M* k (B = 0) = ~ Jo , (34) 

where we assumed the symmetry that the left and right 
coupling are equal to the fixed but arbitrary value Jo- 
Note that this model has been studied in detail by one of 
the authors in Ref. [13 where the system showed a current 
through the left dot even when voltage was applied on the 
right side (transconductance effect) as discussed in more 
detail at the end of this work. 



III. FLOW EQUATIONS IN LOWEST ORDER 
A. The method of flow equations 

In a model with a clear separation of energy scales it 
is often an advantage to focus on the low-energy physics 
and find an effective representation for the high-energy 
physics. In the poor man's scaling approach as will also 
be discussed in more detail in the next section, the scat- 
tering processes containing energies at the large band 
edges are integrated out. Their effective contribution is 
put into a renormalized interaction and by further reduc- 
ing the band cutoff one arrives at an effective model for 
a low-energy band. 

The separation of energy scales is also important in the 
flow equation method. The philosophy here is again to 
find an effective Hamiltonian which describes the same 
physics as the original model but can be easily solved 
like a quadratic Hamiltonian. As an example the kinetic 
Hamiltonian e/so-ct a c ka is diagonal in the conduction elec- 
tron momenta, whereas the interaction part proportional 
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to : Sk'k — a iT a ' crC-ko connects electrons with different where 
momenta. In the matrix representation of the Hamilto- 
nian H we separate the diagonal contributions, H , from 
the off-diagonal contributions Hi nt . The aim of the flow 
equation method is to generate an effective H(B) start- 
ing with some general flow parameter B = and modify 
H{B) accordingly such that H(B = oo) is diagonal. 

In the flow equation method we achieve this by in- 
finitesimal unitary transformations 21 



dB 



H(B)=[ V (B),H(B)]. 



(35) 



= -Tit 



Note that this expression forces the generator rj 
to be anti-hermitian which is equivalent to claiming that 
the transformation is unitary. The choice of the gener- 
ator can be different from problem to problem, but the 
canonical generator- 



's) = [H ,H int (B)} 



(36) 



has proven to be a stable choice. The canonical generator 
automatically fulfills = —77. Since it is the product of 
two Hamiltonians it is proportional to energy 2 . Conse- 
quently the flow parameter B is of the order of 1 / energy 2 . 

In contrast to the standard scaling theories this choice 
of rescaling eliminates the outermost components in the 
matrix Hamiltonian or in other words: scattering pro- 
cesses involving an energy transfer of the order of the 
frequency band-cutoff A are integrated out in the course 
of the flow procedure, B = 1/A 2 . The final Hamiltonian 
contains only energy-diagonal processes with a renormal- 
ized energy. It is a clear advantage to keep information on 
all energy scales, in particular for the non-equilibrium sit- 
uation where scattering processes away from the ground 
state play an important role. 

For a more extended introduction to the flow equation 
method we refer to Rcf. m As a side remark we want to 
mention that corrections from taking into account normal 
ordering with respect to the interacting ground state are 
of fourth order in the interaction^ and can thus safely be 
neglected in our calculation to third order. A first appli- 
cation of the flow equation method to coupled quantum 
dots in equilibrium can be found in Ref. |22L 



B. Flow equation for the double quantum dot 

system 

For the double quantum dot system studied in this 
paper the generator r\ is chosen to be the canonical gen- 
erator 77 = [ifoi-fJint] and given explicitly by 

V = E XXT J : + SR)s(k'Mkj) : ( 37 ) 

j=L,R fc'fc 

+ E E^ MJ = (0L - S R ) ± 2i{S L x S R ))s [k , j){kj) 

j=L,R k'k 



sum. j 
Vk'k ' 

Pj 

Mj 
Vk'k 



tk > J k'k > 

e k -K)Ml k . 



(38) 
(39) 
(40) 



As discussed before we can here observe that the coupling 
P k , k or M 3 k , k corresponds to a transition between a singlet 
and triplet state with an energy cost of ±K, respectively. 

Due to the construction of the canonical generator (oc 
energy 2 ) the flow parameter B is related to the tradi- 
tional energy/frequency cutoff A by 



B cx 



A 2 ' 



(41) 



Inserting the canonical generator 77 into the flow equa- 
tion, Eq. (|35p . we find in lowest, linear order an expo- 
nential behavior of the coupling functions. Thus we can 

define an effective coupling J k r k ,:> (B) 



J s k T 3 (B)=e- B ^'-^ JlJr(B) 



(42) 



and 



Pi k (B)=e- B ^'-«* + V pi ik{B ), (43) 
ML(B) =PUB) =e- B ^'-<>-V 3 pj£(B), (44) 



The effective couplings J™™' 5 and P k , k obey a scaling 
equation with a scaling function /3 which has to be de- 
termined from higher than linear order terms. The expo- 
nential dependence though mirrors the physical picture of 
the Kondo coupling: It is logarithmically divergent when 
energy scattering processes with initial state k and final 
state k' are energy-degenerate, e.g. e~^ ek '~ €k ^ A " > = 1 
for ey = tk , and away from the coherence conditions the 
coupling functions are suppressed, see also Fig.[2]and cor- 
responding discussion. In contrast to the single-impurity 
Kondo model (without magnetic field), a divergent cou- 
pling for P k , k representing the singlet-triplet transition 
can only be expected when a scattering process in the 
leads matches the energy of a transition inside the quan- 
tum dot. 

The full expression for the flow equation calculation 
to second order in the interaction, the so-called one-loop 
order, is given in the appendix. In Fig. [5] we show full nu- 
merical calculations for these one-loop result in the case 
of k' = k and compare with the solution obtained by 
the diagonal parametrization^iS,. This is a by now well- 
established approximation that allows some analytic in- 
sight into the flow equations and simplifies the numerical 
effort significantly. 

In the diagonal parametrization we assume that the 
important energy dependence €k (momentum k) is given 
by the exponential decay e~ B ^ k '~ tk+aK ^ and we can 
approximate 



-B(e k , 



- tk+aK)2 f(ek',e k )^e- 



-B(e k ,-e k +aKy 



(45) 
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where a = {0, ±1} and = (e k ' + £fe)/2. 

Starting from two energy arguments for the incoming 
and outgoing conduction electron one energy is kept fixed 
but arbitrary and the other is assumed to fulfill the equa- 
tion ey — e k + otK = 0. For example for the coupling to 
the total spin (Sl + Sr) this yields: 

J^r(B)=e- B ^-^J™r4 y2 (B) (46) 
where J^(B):=^(B). (47) 

In the diagonal parametrization for Pli k /Mli k one has 
to be cautious since the choice ey — e k = —K in P k , k is 
not unique. The correct momentum dependence is only 
recovered if we choose the diagonal parametrization as 

PUB) = e- B ^-^+*) 2 P( k , +k)/2 {B). (48) 
where P* k {B) := P J ek _ K/2>ek+K/2 (B) (49) 

Note that in the definition of P/. the average energy 
is given by e s = (e k - K/2 + e fe + K/2)/2 = e k . For 
the coupling M 3 k , k = P kk , the assumption in Eq. (1451) . 
i.e. ey — K/2 = ek + K/2, is automatically fulfilled and 
in diagonal parametrization: 



L k'k 

since 



(k'+k)/2 



(B) 



(50) 



Pl{B)=Ml k+K/2tk _ K/2 (B) = P- t 



k -K/2,e k +K/2( B )' 
(51) 



Using the assumption in Eq. (|45j) which leads to the 
diagonal parametrization we arrive at the one-loop flow 
equations 



dB 

-Y J {l~2n{ qj ))(e k -e q )e' 2B ^-^\j\ 



sum j \ 
(k+q)/2> 



4^(1 - n(qj))(e k - e q + K)e -^ k -e q+ Kf ^ 



k+q)/2l 



(k+q)/2) 



(52) 



and 

dPj{B) 
dB 



- - n(qj)) {2{e k -e q - K/2)) e -™^-«-*/*)' 

q 

J {k-K/2+q)/2 r (q+k+K/2)/2 

+ ]T n(qj) (2(e k - e q + K/2)) e -^-^+^f 



pj 7-sumj 
r (k~K/2+q)/2 J (q+k+K/2)/2 



(53) 



For details of the calculation we refer to the appendix El 
Note that M{, k is given immediately by the solution for 
P k in the diagonal parametrization. 

The one-loop order contains the integration over one 
internal degree of freedom, the momentum q. Assum- 
ing a constant density of states (DOS), N(0) = 1/(2A ), 
of a fiat band with bandwidth Ao around the Fermi en- 
ergy, we transform the summation over momenta q to 
an integral over the energy e q : J2 q ^(0) J*_a ^e q - 
The DOS is absorbed into the dimensionless couplings 
9k ^ = N(Q)Jl^ ^d p> k , k = iV(0)4 fc . 

Predicting that the couplings will be only logarithmi- 
cally dependent on the energy we assume that the energy 
dependence is dominated by the exponential function and 
simplify 

f(x) exp(-2B(a; - cf) /(c) exp(-2B(a; - cf) (54) 



Then the integration in Eqs. (J52J) and (|53[) is indepen- 
dent of the energy argument of the coupling function. 
Straightforwardly we find for example 



de g (l - 2n(qj))(e k - e q + aK)e 



'j e -2B( tk -V J /2+aKy 



1 + r,- 



-2B(e k + V j /2+aKY 



(55) 



where we used the two-step Fermi function, Eq. (|29[) , and 
a as placeholder for the corresponding prefactor to K. 
If we assume that the leads are symmetrically coupled 
(asymmetry parameter r,j = 1) we can write the flow 
equations as 



dg s r^{B) _ 1 



1 



dB 



y - e 

25 ^ 2 

v=±l 
n 1 V- 1 

+ 2 — > - 

2B ^ 2 

u=±l 

+ 2 2B E 2 

u=±\ 



2B(e k +^V J /2) 2 fsum,j\ 

-2B( tk +K+uV J /2) 2 ( j V 
1 \t'k+K/2 ) 

-2B(e k ~K+uV 3 /2) 2 ( j V 
e [Pk-K/2 J 

(56) 



and 

dpj(B) 
dB 



J_ y l_ 
2B t-{ 2 

— V - 

^ 2B 2-f 2 

v=±l 



p -2B(e k -K/2+uVj/2Y sumj j 
e 9k-K/lPk 



-2B(e k +K/2+vV J /2Y J sumj 
e Pk 9k+K/2 



(57) 



Note that here the term at e k = is exponentially small 
for B > (Vj/2) 2 due to e~ B ^/ 2 ) 2 (assuming K = 0). 
Reducing the band cutoff A to in a system with a large 
applied voltage thus leads to a failure of the theory. If 
the voltage is taken beyond the linear response regime 
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0.25 




0.15 



FIG. 2: Flow of the coupling <?^™' L for ty = in the 
numerical exact calculation (dashed lines) and for g k uln in 
the diagonal parametrization (solid lines) versus energy e& 
for three different values B = I/A 2 and A = K/3, K, 3K. 
Further parameters of this plot are chosen K/2K.Q = 0.05, 
and gr(Ao) = 0.128 leading to T K = 0.01. Note that the flow 
of K is neglected here. 



it is necessary to study the frequency/energy-dependent 
behavior of the coupling functions and the divergence at 
the two new Fermi edges e& = ±Vj-/2. 



C. Discussion of results 

Note that the left and right leads do not mix in 
Eqs. (|56]) and |57|) and therefore the left and right cou- 
pling can be studied individually. 

In Fig. [J] the flow of the coupling g s k ™' L is plotted 
for ek = tfc' versus the energy ek for different values of 
the flow parameter B while the voltage V is set to zero. 
Since there is no cross-coupling the scaling behavior of 

sum.ii! • i i ■ i , sum.L mi • i 1 r sum.I 

9k' k 1S identical to g k , k . Ine initial value ot g k , k 
is energy-independent but very soon a frequency depen- 
dence is produced, and for large enough B peaks at the 
Fermi edge efc = and at non-zero energy = ±K 
are visible. Values away from = 0, ±K are exponen- 
tially suppressed. As obvious from Fig. the diagonal 
parametrization g s k um ' L reproduces the frequency depen- 
dence of g s k ^' L perfectly and also the absolute value of 
the coupling is in very good agreement. Thus we can 
conclude from the numerical comparison in Fig. [5] that 
the diagonal parametrization is a good approximation. 

This has proven to be the case for most problems 
studied previously in the context of non-equilibrium 
Kondo models solved with the flow equation method^. 
In the following we will also show that the diagonal 
parametrization is equivalent to the non-equilibrium scal- 
ing method by A. Rosch et ali£ and can thus also confirm 
the accuracy and the validity of the assumptions in the 
other approach. 

For B <C min[l/if 2 , 1/e 2 , 1/V^ 2 ] the scaling equation 
for gj™'"' and p k are identical since the exponential 



e -2B(e k + a K+vV 3 /2f _ 1 does nQt cutoff the flow . Note 

though, that the initial value o{p> k (B ) = 1/2 g s k m ' J (B ) . 
Using B = 1/ A 2 where A is the frequency cutoff the flow 
equation (|5S|) thus simplifies to the well-known Kondo 
coupling scaling function for a spin-1/2, 



dg 
din A 



(58) 



This equation gives us a one-loop Kondo temperature of 
T K = A c- 1 / 29 ^ . (59) 



Thus we expect the flow equation (|58|) to diverge and 
the lowest order to break down when B reaches the value 
Nonzero temperature and current induced deco- 
herence can remove this divergence as we will explain 
later on. However, first we want to introduce the pertur- 
bative RG approach in the next section. 



IV. GENERALIZED PERTURB ATIVE RG 

The interaction of a Kondo problem has the general 
structure 



na,ma' 
1,1 



(60) 



where 7,7' refer to the eigenstates of the eigenstates of 
the double quantum dot (singlet and triplet in this ex- 
ample), m,n are the lead indices, a, a' are spin up and 
down states and T a i a refers to the Pauli spin matrix and 
Ty 7 is a generalized Pauli matrix which has to be chosen 
in accordance with the internal structure of the quantum 
dot setup. The momentum dependence of the interaction 
is neglected and thus we write the Kondo interaction in 
terms of the momentum integrated conduction electron 
operators, c na = Ylk c nka- 

This Hamiltonian can be furthermore described by a 
general vertex 



V. 



tict,u c ;mcr ,w 



7to\u? c :m<7 ,uj „ 



(61) 



Note that the Hamiltonian necessarily has to be hermi- 
tian and thus 



if all couplings are real numbers. 



(62) 



A. General scaling equation 

The idea of scaling has already been introduced. In- 
stead of calculating a perturbative series of a physical 
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quantity, we do a perturbation in the couplings of the in- 
teraction. This provides us with results beyond standard 
perturbation theory, which is known to fail in the Kondo 
model. 

One of the first scaling theories was introduced as 
"poor man's" scaling by P. W. Anderson in Ref. [23| and 
generalized to non-equilibrium by A. Rosch et ali^. Im- 
portantly, the physics of the problem at hand, i.e. the 
expectation value of any physical observable, needs to be 
invariant under scaling. Anderson studied the T-matrix 
for the Kondo model and showed that a scattering into 
a high-energy state can be absorbed into a lower-energy 
setup by including the process to the bandedge A and 
back into an effective interaction for a smaller bandwidth 
A — dA. There are two contributions equivalent to a scat- 
tering process to the upper band edge (electron = Cooper 
contribution) and lower band edge (hole = Peierls con- 
tribution)^. By successively reducing the band cutoff A 
in infinitesimal steps dA we generate an effective inter- 
action at low energy scales which is of the same form as 
the original one but contains a renormalized Kondo cou- 
pling J (A). The change is described by a scaling equa- 
tion dJ(A)/dA and including only one virtual state at 
the band edge is referred to as one-loop order. As the 
poor man's scaling is also perturbative, a truncation has 
to be made and only certain renormalization diagrams 
are taken into account. 

In the traditional poor man's scaling^ the band cut- 
off is reduced to zero, A — > 0. This is a problem in a 
non-equilibrium situation as emphasized before because 
energies beyond the ground state play an important role. 
In Ref. the approach from Anderson was generalized 
to renormalize only one of the two band-cutoffs, e.g. the 
outgoing energy hu>' c for a vertex while the incoming en- 



ergy Hw c is fixed (but arbitrary). In this context the 
Cooper and Peierls contribution to the scaling equation 
have to be calculated in Keldysh notation and we find 15 
that the leading logarithmic contribution originates from 
an integral of the form 

A 

tJ-t I de-^—ta,nh(-^=) « ^— (63) 

<91nA J x-e \2TJ x-A v ' 

-A 

i.e. the real part of the quantum dot Green's function 
and the lesser part of the conduction electron Green's 
function including the sharp step of the Fermi function 
at zero temperature. The external energy x can be the 
exchange energy K, magnetic field h, chemical potential 
fi or combinations of those. 

In the perturbative RG method we approximate 
A/(A - i) by 0(A - |a;|) since A/(A - x) « 1 for A > x 
and A/(A - x) « -A/x for x > A. At this point 
we like to state that one major difference between the 
two scaling method is the choice of the cutoff-function. 
In the flow equation method the cutoff function e~ Bx 
is valid in general and therefore we expect a better res- 
olution in the proximity of the logarithmically enhanced 
peaks. For further information on the pRG method we 
refer the interested reader to Ref. |^ or [lfl 

In the general notation 

= £ E K"^' d A^ (64) 

77' ncr^mcr' 

one can derive a scaling equation of the form 



,ncr,cd c ;m<7 ,aj„ 



A=±l r],v,s 



where we introduced the notation = 0(A — |x|). 
Eq. (I65p is the generalization of the equations given in 
Ref. 0- A derivation can be found in Ref. [TH Every in- 
going and outcoming leg of the vertex is assigned with a 
frequency, w c ,uj' c and for the conduction electron 

and quantum dot pseudo fermions, respectively. 



B. pRG for the Double Quantum Dot System 

We now concentrate on the case of the double quan- 
tum dot system where the general vertex is of the form 



Eq. (|6II) . Furthermore we assume that there is no ex- 
ternal magnetic field applied and therefore the setup is 
spatially invariant, i.e. the three triplet states t+,to,t- 
are degenerate. Like in the flow equation approach there 
are thus only three vertices: a triplet-triplet transition 
without energy cost/gain in the DQD, a singlet-triplet 
and a triplet-singlet transition involving such a process. 

The general vertex has four frequencies assigned to it, 
where one frequency is fixed due to energy conservation, 
which we impose on the vertex. The quantum dot is de- 
scribed by pseudo-particles d 7 which have to obey a con- 
straint. In the following we set the energies "on-shell", 
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e.g. the frequency of the incoming particle w 7 is given by 
the eigenenergy of the state 7, i.e. e 7 neglecting a finite 
lifetime due to hybridization with the leads: 



V. 



na,cj c ;m(7 ,u c 



V 



(66) 



With the energy conservation the vertex thus only de- 
pends on one frequency which is chosen to be the ingoing 
frequency in the following 



V 



ncr.u} r ;m<j .uj^ 



7,e T ;7',e 



77 < K). 



(67) 



Evaluating now the product of Pauli matrices we ar- 
rive at the three scaling equations for the dimensionless 
Kondo couplings g™y (w) = A r (0)J"™(w) for zero mag- 
netic field 



din A 
<91nA 
<91nA 



2 ^ 



(2 S S"(o;)^( W )e ( , 



+2 ff Sf (w)^(w + JQO^+k) , (68) 

(2 S ™( w -if)3 s TMe a -„- K 



1 

2 ^ 

V 

+2g%(u))gT(u)& n ), (69) 



+2^(w)^(w)e < ,- / ,„). 



(70) 



See reference [l5j for the details of the derivation. 

Note that in order to arrive at this scaling equations 
the following approximations had to be made. First, the 
cutoff is sent to on the right hand side of Eq. ([65]). 
Otherwise the imposed energy conservation assumed for 
the left hand side is not fulfilled on the right hand side 
and the RG equation is not self-consistent. Second, the 
pseudo-fermions describing the quantum dots states are 
assumed to be on-shell, Eq. (pjB")) . This step is not nec- 
essary in the flow equation approach since the impurity 
spin is kept as an operator without introducing pseudo 
particles. Third, the energy conservation on the vertex 
is enforced and thus we end up with only one energy in- 
dex (which is chosen to be the ingoing energy). Fourth, 
the approximations in Eqs. (|45p and ([51)) are used in the 
integration identical to the approximations used in the 
flow equation method. 

Note that the hermiticity of the Hamiltonian which 
leads to pk'k — mkk' corresponds to 



9 st M = 9ts 



(71) 



in the perturbative scaling approach. 

In order to be able to compare the two results we in- 
troduce a new symmetrized coupling g"™(cj) analogous 
to the diagonal parametrization of p J k which is defined as 



nm { 
ht 



» = 2£T( W - K/2) , 
£» = 2£> + X/2), 



(72) 
(73) 



which thus fulfills the same initial condition as p 3 k 

££/>)Ia=Ao = i<#> + K/2)\ A = Ao = IjV(O) Jo (74) 

and is peaked at ±K/2 instead of at 0, K and —K, 
where g"™ and g"™ show resonant features. 

With this new definition and inserting the approxima- 
tions as mentioned above yields the scaling equations 



<91nA 



u-Hv-K/2 



+g™(u - K/2)~gT(")®u 



(75) 



<91nA 



£ (2#J> + K/2)g%(u, + K/2)Q n+K 



+2&> - K/2)gT(u> - K/2)G U ^- K 
+g™{w)g%(u)Q n ) (76) 

In most applications a further convenient approximation 
is used: in the frequency integral over all coupling func- 
tions it is assumed that the main contribution arises from 
the value at which the cutoff functions vanish. Thus the 
set of equations reduces to a parametric set instead of a 
continuous function. 

V 

+9?t V M9T(»u - K/2)Q u _^ +K/2 ) , (77) 



9 In A 



PaTOh* - K/2)~g™{^ - K/2)e n+K 

V 

+2g%( Jh , + K/2)g»?(»v + K/2)Q U _^_ K 

+gT{»X(^)®*-^) (78) 



We leave out this step while comparing the expressions 
([75]) and {75]) directly with Eqs. {55]) and {57]) derived 
with the flow equation method. However, the latter ap- 
proximation is used in the numerical routines to acceler- 
ate the calculation. 



COMPARISON OF THE TWO METHODS 
TO ONE-LOOP ORDER 



In the following we will show that the two methods use 
the same approximations and therefore are identical to 
leading logarithmic order. 

It is straightforward to see that the couplings in the 
two different calculations are related by gtt{w) ~ <?| um 
and g s t(w) ~ Pk, where the energy ek is to be identified 
with the frequency (energy) ui (huj). 

We rewrite Eq. ([55]) and Eq. ([To]) to study in detail 
the commons and differences of the pRG and the flow 



10 



equation method 



E 



+ 2 + K/2)gX(u + K/2)@ u+K .^ 



K/2)g™(«i - K/2)Q U . K . 



dln(B-V2) " 



2 S 



sum, 7 

0* 



1 -2B(e t - |U „) 



2 (^+K/ 2 ) r 
(Pfc-JC/2) 



-2B{t k +K-^ v ) 2 



-2B(e h -K-n„ 



(80) 



2 ^_k/ 2/ 

The flow parameter _B is related to the frequency cutoff 
by £? = A~ 2 and fi v — ±Vj/2. Therefore the equations 
are actually identical down to the prefactors. The pref- 
actor 1/2 in Eq. (ISTjf stems from the even-odd combi- 
nation of the leads where g even = (1 + r)g nm and thus 

g even = 2g nm f or r = 1. 

Let us repeat the common approximation before we 
embark on the differences. Both methods take into ac- 
count only the leading logarithmic order by treating only 
the lowest order diagrams which contribute to the scal- 
ing. In the perturbative RG the pseudoparticle energies 
are taken to be onshell which corresponds to treating 
the spin operators without a bath. In both methods the 
level of complexity is reduced by imposing energy conser- 
vation on the vertex, i.e. the energy of the outgoing con- 
duction electron is given by the energy of the incoming 
electron diminished by eventual inelastic processes inside 
the dot. Since both methods keep information on the 
whole bandwidth and the renormalization reduces only 
the bandwidth of the state that one scatters into, they 
are destined to treat non-equilibrium physics on different 
energy scales, i.e. when physics of more than the ground 
state play a role. 

The main difference are the two different cutoff func- 
tions e- Bx2 = e - {x /^ 2 versus Q x = 0(A - The 
step function is an approximation in the perturbative RG 
whereas the exponential cutoff arises naturally in the flow 
equation method approach. The form of the cutoff influ- 
ences at most the lineshape in the proximity of logarith- 
mically divergent couplings, see Fig. [3] and 01 Since these 
are heavily influenced by higher-order renormalization, 
their shape in leading logarithmic order is not reliable 
anyway and only the flow at the corresponding energies 
is described correctly. Note though, that for example in 
the current or other physical observables, it is the aver- 
age over a voltage window which enters, i.e. the integral 
over the frequency-dependent couplings. Therefore this 
difference can play a role when comparing result for the 
non-equilibrium current. 

Furthermore there is a summation over the lead in- 
dices in the perturbative RG. In the flow equation calcu- 
lation we used the symmetry that only the even chan- 
nel of each lead is coupled to the dot and therefore 



the lead index drops out. In general the model can 
also be studied without doing the even-odd transfor- 
mation 7 , e.g. to study models which are not equivalent 
to the Anderson impurity model by the Schrieffer- Wolff 
transformation. As explained in the introductory part, 
e.g. the important transport part g 12 between two leads 
on the left side is given by the asymmetry parameter and 



(79) g L ' R ' L = (l+r L )g L / R < 22 : g L / R < 12 = y/rE/0-+r L )g 



L/R,L 

We concentrate for the comparison only on the symmet- 
ric case r L = 1 and consequently find that gl um = 2g tt (u>) 
and p k = 2g st (u)), 

The strength of the flow equation method is that the 
diagonal parametrization is a convenient choice in order 
to find analytical expressions, but not a necessary limi- 
tation of the method. In practice the numerical cost lim- 
its the evaluation of the flow equations to the diagonal 
parametrization. 

Nevertheless both methods are still bound to break 
down due to the strong-coupling behavior of the Kondo 
correlations. In the following we will explain how deco- 
herence effects are included in the two different methods, 
how they provide an additional cutoff and in the end 
compare the two methods again. 



VI. BEYOND ONE-LOOP 
A. In perturbative RG 

The Kondo problem at zero temperature is known not 
to be solvable by perturbation theory. Even though the 
renormalization schemes can improve the limit of validity 
due to an appropriate summation of diagrams, the theory 
still breaks down at the energy scale A = Tr-. However, 
in equilibrium at finite temperature T > Tk or in non- 
equilibrium for sufficiently large voltage bias one expects 
a well-behaved weak coupling expansion. 

The seminal work of A. Rosch et al. 2 on the non- 
equilibrium Kondo model achieves this by including a 
physically motivated cutoff given by the current-induced 
noise in the system. Decoherence is unavoidably present 
due to the non-zero steady current flow through the sys- 
tem, which generates Johnson-Nyquist current noise. In 
other words, the quantum dot states gain a finite lifetime 
due to elastic or inelastic cotunneling processes with the 
leads. Also decay rates due to external baths can play 
a role. In general a non-zero decoherence rate Y has to 
be included in the retarded Green's function. This leads 
to a new scaling behavior since the leading logarithmic 
diagrams changes as 



a£- I ,h 



x — e , , / e \ 



■ tanh ( - 



dA J ~~{x-cf+Y 2 \2TJ 

-A 

x - A , x + A 

"A-; r~To ^77 + A- 



(x-A) 2 + r 2 (.T + A) 2 + r 2 ' 



(81) 



(82) 



11 



Thus the cutoff has to be corrected to be 

e x = e(A- Vx 2 + r 2 ) (83) 

instead of just 0(A — \x\). It was shown perturbatively 
in Ref. [lH that in non-equilibrium selfenergy and vertex 
correction become important to the same degree. Fur- 
ther studies using different rcnormalization group meth- 
ods confirmed the observation that the non-equilibrium 
dccoherence rate are determined by transport processes 
and can be different to the thermodynamically expected 

4 5 24— ?fi 

expressions i i i — 

The statement that decoherence terms have to be in- 
cluded is equivalent to the failure of the on-shell assump- 
tion since the levels in the quantum dot are broadened 
and thus the spin state on the quantum dot gains a finite 
lifetime. To find the correct cutoff one has to calculate 
the spin susceptibility and find the correct Lorentzian 
shape. 

In Fig. [3]we show the flow of the triplet-triplet coupling 
<7™ t m (w) for n = m = 1 as a function of the frequency w 
and for different values of the cutoff A. 



produces this behavior automatically, which is one of the 
advantages of the pRG method. 

Without the cutoff given by T the couplings at the 
Fermi energies uj = ±V/2 are logarithmically divergent 
and the poor man's scaling approach breaks down. In- 
cluding r in the cutoff function 0r simulates the phys- 
ical observation that similar to a non-zero temperature 
T ^ the infinite series of infinitesimal spin-flip excita- 
tions around the Fermi edge is interrupted by a non- 
coherent process. The flow of the coupling stops at 
gtt oc l/ln(r/Tff)^ assuming that F is the relevant cut- 
off (which can otherwise also be proportional to T or 
combinations of internal energies). 

Note that the divergence of the coupling at the Fermi 
edge is not necessarily equivalent to the breakdown of 
the calculation. For example in the calculation of the 
non-equilibrium current the average of the frequency- 
dependent coupling over the whole voltage window needs 
to be included where the divergent coupling is just a 
boundary term2&. 



B. In the flow equation method 




A=2V 



025 1 1 1 1 1 1 1 1 1 

TB.2 -0.1 0.1 0.2 

CO 

FIG. 3: Poor man's scaling including a frequency de- 
pendence of the coupling for a Double Quantum Dot 
system; gU = 0.025, N(0)K = 0.025, T « 0, 
N(0)(eV L ) = V = 0.1, N(0)(eV R ) = 0; values of the cutoff 
A/V = 2, 1, 1/2, 1/5, 1/10, 1/15, 1/25, 1/50, 1/100 where the 
first (solid line) and the last value (dashed line) are marked 
extra in the plot. 

For a large cutoff A the coupling is a constant and does 
not depend on the frequency of the scattered electrons. 
As the cutoff is reduced a frequency dependence evolves 
and the coupling continues to grow at the resonant energy 
scales, i.e. at uj — ±V/2, ±K ± V/2, while the flow does 
not continue at other energy scales due to Q w . Once the 
energy scale which is needed for a scattering process is 
smaller than the bandwidth, no further processes can be 
integrated out and it is expected that the flow should not 
continue. Including the frequency dependence and the 
frequency cutoff W in the poor man's scaling equations 



In the flow equation method, the generalization to 
higher orders in the coupling function is straightforward 
as has already been illustrated in Refs. 

In second order in the Kondo coupling J new interac- 
tion terms are created. These can be compared to a 6-leg 
vertices with two incoming and two outgoing conduction 
electrons and one incoming and one outgoing spin state of 
the quantum dot in the pseudo-fermion language. Con- 
tributions to the flow of the four-leg Kondo interaction 
vertices are generated by a product of one of the newly 
generated vertices with the initial Kondo vertex. Note 
that since two conduction electron lines are thus inte- 
grated out, these terms are very similar to the self energy 
and vertex corrections discussed in the previous section. 

In contrast to the one-loop contribution, the two-loop 
contribution includes an integration over two interme- 
diate energies, i.e. an electron-hole pair like n(qv)(l — 
n(q'v)). This integration is taken over a window of 
energies, e.g. voltage window or singlet-triplet transi- 
tion window. Besides the standard two-loop term in 
the /3-function dg/dlnA = —2g 2 + 2g 3 for a spin-1 we 
also expect to find contributions from g 3 V/A or g 3 K/A, 
i.e. coupling 3 x phase space. These arise naturally in the 
calculation and can be identified as the cutoff rates given 
by e.g. the Johnson-Nyquist current noise. 



1. Discussion of the flow to second order 

In the case of a double quantum dot system these newly 
generated couplings are 
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H int = ^2 K k^k-q'q i : (Sl + Sr) (s(k>j)(kj) x «(?'iO(gto) : 

k'k;q'q 

+ ^2^2 K k'k;q'q i : ({Sl - Sr) + 2i(S L X Sr) J (s(k'j)(kj) X S,y „)(„,) 
ji; k'k\q'q 

+ XI X! K k'k-q'q i ((Sl ~ Sr) - 2z(S L X S H )J [S(k'j)(kj) X S(g'v)(fftO 
iu k'k;q'q 

+ J! H K £k *q i2i ( : (-^L X S(k'j)(kj)){SR X S(g'«)(5t,)) + (-Sfl X S(k'j)(kj)){SL X %'*)(<po' 
ju k'k;q'q 



(84) 



r 



These interactions are not present in the initial flow 



K ! Wk\q'q\B 



0) = 
= 



= 0) 

= 0) = 
K^{B = 0) = 



(85) 
(86) 
(87) 



They can be interpreted as an entanglement of the quan- 
tum dot states with the spin states of the conduction 
electrons. This can be seen by studying of the flow of 
the coupling function where we find that the initial free 
spin state (Sl + Sr) evolves into the entangled spin state 

(Sl + Sr) x s/q' v )(qv) f° r B — > oo. Correspondingly we 
observe in the flow of the Hamiltonian that the coupling 
Jk'k starts to decrease at an energy scale when Kk>k;q'q 
starts to grow. The flow of Jk'k is reversed in the sense 
that it does not diverge logarithmically but for B — > oo 
we observe Jk'k — ► while the values of the couplings in 

H int S row - 

Note that also an RKKY-like interaction is created to 
this order. All kinds of potential scattering contribu- 
tions are neglected here, similar to the pRG calculation, 
because they do not contribute in the wide band limit. 
Since we study only antiferromagnetic coupling between 
the two spins we can neglect the RKKY interaction in 
the following. For a discussion of the flow of this cou- 
pling we refer the interested reader to appendix [Bl where 
it is shown that the coupling K is similar to a magnetic 
field in a spin- 1/2 Kondo problem renormalized by the 
coupling to the leads^. Note that the change of K is in- 
cluded in the numerical results shown here if not stated 
otherwise. 

The rather complicated entanglement term 

K WkrfJ i2i ( : (Sl X S(fc'j)(fci))(^R x &(q'v)(qv)) 
+ {Sr X %j)(fej)) {S L X S { q'v)(qv)} 
= (-4)K$f. q ,q : (SlSr) {s(k'j)(kj)S(q'v)(qv)) 



2KS 



■;q'q 
RKKY 



(SLS {k ' ] )(k ] ))(SRS {q , v){qv) ) 
2K k'k-q'q '■ (SRS(k'j)(kj)) {SLS(q>v)(qv)) 



k'k\q' q 
-RKKY 



(89) 



reproduces an RKKY-like interaction. Note that during 
the flow S L S R -> (S L x s^^^Sr x s [q , v) ( qv) ) and 
therefore the additional term is effectively a rescaled ex- 
change interaction. This term is not present for example 
in the spin- 1/2 model and special to the exchange cou- 
pled quantum dot. 



We now study systematically how the new terms in the 
Hamiltonian decompose into the initial coupling terms. 
A detailed expression can be found in the appendix [C] 
Here we write down only the final expression for the sym- 
metric setup discussed throughout the paper. The flow of 
the coupling g^™' 3 (B) in diagonal parametrization yields 



dgT m ' 3 (B) _ 1 1 



dB 



9R 9 



-2B(e k +vVj/2) 2 ( sumj 



2B2 
1 1 



/ sum, i\ * 

\9k ) 



v=± 



+ 9R 9 2^ e \Pk+K/2) 



v=± 



+ i~Y,e- 2B(ek+vV > /2 - K)2 (pi-K / 2) 



v=± 



-^Y.9 S r d (B){9Zaf{B)f 

V 

-^W,V v )g k ™' j (B)(g°™f(B)) 2 

V 

_ (J_ C -2BK> k ^ , 

K 4B 4 V 2B v " 

xJ29 s r d (B)(2p v max (B)) 2 

V 

-Y,f{K,V v ) g^(B)(2p° max (B)f , 



(90) 
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and 



by 



dpi(B) _ 1 1 



dB 2B 2 



v=± 



-2B{t k +vV ] /2~K/2Y sum j i 
yk~K/2l } k 



1 1 



_iV -2S(£ fc +fV,/2+_ftr/2) ;i T) J„sumj 
9R 9 



-^/(o,KK(B)(C"(fi)) ; 



1 „-2BK 2 K I TT 



, —J — erf(V2Bif)) 
x^p{.(B)(2^ a:c (B)) 2 

^/(ir,y„y fe (s)(2^ a:c ( J B)) 2 ) (9i) 



where for = tr = \ 

= ^^(e- 2B(X " K)2 +c- 2B ^ +v ^ 2 - 2c- 2BK2 
+ \\[^\ [i K - V v )erf(V2B(K - V v )) 



+ {K + V v )eif{V2B(K + V V )) 
- 2Kerf[V2BK)] , 



(92) 



and in particular 



/(o,v.) = -Li(.-* 



i 



(93) 



Note that we assumed that the frequency dependence 
is given dominantly by the exponential decay and there- 
fore the couplings in the integration over the energy 
window can be approximated by their most divergent 
term, i.e. S(™'* )/2 ~ 9^ v = 9^ aa ' v if voltage V v 
is zero. For V v = there is no contribution from 
/(0,0) = f(K,0) = 0. On the contrary for finite volt- 



age and B > V v 2 : /(0, K) ~ ^ 2 V 2B ' 2 V 2B 
since 3> 55 an d /(0j provides the leading con- 

tribution to the flow. 



2. Two-loop results of the Double Quantum Dot system 

In the equilibrium case, where the applied voltages 
Vl = Vr = 0, the scaling equation for g^ xm is given 



dg a r J (B) 
dB 



-2Bei ( sumj 



2B 
1 



\9k ) 



2 — c 
2B 

1 

4B 



-2B(e k +K) 



(Pk+K/2) 



-2B(e k -K) 



Pk-K/2) 

\Y.9 S r' i {B){9 S rZf{B)f 



[ie- 2Ml + -*f^eri(V2BK)) 
K AB 4 V 25 v " 

x^ 5 ™'( J B)(2^ ox (S)) 2 . (94) 



Note that the scaling equation is of the form of a two- 
loop calculation with an additional contribution pro- 
portional to (K/ V 2B)ei~i(V2BK) . For large values of 
the argument in the error function we can approximate 
erf (2) — > sign(a;). Therefore we find in the limit B — > 00 
that this term is oc 1/ \/B > l/B and thus this term 
dominates the asymptotic flow. As discussed later on, 
this provides a cutoff in the flow given by the singlet- 
triplet exchange interaction K. This cutoff can prevent 
the divergence of the elastic Kondo coupling g| um by fa- 
voring a singlet formation of the two dots instead of a 
Kondo singlet with the leads if K is larger than a critical 
value. This is a signature of the quantum phase transi- 
tion in our system and will be discussed in more detail 
in another publication^. 

Furthermore we try to illustrate the effect of these de- 
coherence terms in the two-loop calculation in the limit 
Vj S> K and refer the reader to the explicit discussions 
in the literature to non-equilibrium flow equations for 
further reading^. 

Note that we have to distinguish between the voltage 
applied on the left side or right side since there is a sum- 
mation over the left and right lead index in third order 
in the coupling. For now we will assume Vr — and 
only Vl = V 7^ 0. In the current this can lead to the 
interesting effect of a transconductance^ as discussed at 
the end of this section. 

For 2BV 2 > 1 > 2BK 2 such that we can set 
evi{y/2BV) ~ 1 and neglect contributions from er 2BK 
we find the scaling equation 



dB 



dg k '\B) _ J_ 1 -2B(e k +aV/2) 2 / sura,i> 

~ 2B L> 2 [9k ' 

a=±l 



9r°' j (B)(9™% v (B)y 

g s r^BM iax (B)) 2 

(95) 




v=L,R 



In this limit the voltage V dominates the flow for B — >• 00 
and thus provides a cutoff scale for the flow. 
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-0.2 -0.1 0.1 0.2 

FIG. 4: Flow of <?^ um,i as a function of energy eu for different 
values of the flow parameter B — (1/A 2 ) where A chosen 
identical to Fig. [3] Furthermore N(Q)K = 0.025 N(0){eV) = 
0.01, T « and g(Bo) = 0.05. Same parameters as in Fig. [3] 

The actual decoherence rate has to be determined 
studying the spin susceptibility of the system. Neverthe- 
less, it has been shown e.g. for the spin-1/2 Kondo modeli 
that the rate T read off from the correlation function is 
equivalent to the observed cutoff in the scaling equation. 

In Fig. @] we show (similar to Fig. [3]) the flow of the 
coupling <?fc UIn '' 5 as a function of for different values of 
the flow parameter B. As in the pRG calculation the 
coupling is initially equal for all values of the energy and 
a frequency dependence evolves slowly with increasing 
values of B. Note the difference between the two calcu- 
lations: in the flow equation method the two-loop con- 
tribution leads to a decrease in the running couplings as 
soon as the energy scale T is reached. At this stage the 

(2) 

higher order coupling terms in H^ nt start to grow and 
the entangled spin state determines the physics. This is 
seen as a decrease in the initial Kondo coupling g™ m, ' J , 
whereas in the poor man's scaling approach the flow is 
just stopped at this scale and then stays constant. 

It is still an open issue if the coupling exactly at 
B = 1 /r 2 can be used if physical observables are calcu- 
lated, and if the cutoff scheme motivated by a self energy 
cutoff is valid to higher orders. For the parameter set 
chosen in Fig. 2] the peak structure only evolves after the 
absolute value of the coupling goes down. Otherwise the 
overall behavior is very similar to the result using pRG 
as shown in Fig. [3] A direct comparison of results of the 
two methods can be found in the next subsection. 



3. Trans conductance 

A non-zero transconductance is found in the double 
quantum dot system^. If the ground state of the double 
quantum dot system is given by the degenerate triplet 
states (K < 0), we find a zero-bias resonance which 
can be Kondo enhanced. On the other hand if the dou- 
ble quantum dot system is in the non-degenerate singlet 



ground state (K > 0), transport is blocked for voltages 
below a threshold given by the exchange interaction K. 
Once the applied voltage is larger than K, transport is 
also allowed which includes the three triplet states and we 
find an inelastic cotunneling step which is logarithmically 
enhanced by the Kondo correlations. A non-equilibrium 
occupation of the triplet becomes possible due to the 
large applied voltage. For example for a large voltage 
applied on the left side the linear response current on the 
right side is not blocked but there is a non-zero signal. 
This transconductance dlnjdVh was studied in detail by 
one of the author o 14 i 15 . 

The transfer of decoherence from one quantum dot to 
another can also be seen in the flow equation method. 
In contrast to one-loop order where only one virtual pro- 
cess is allowed, there are two intermediate states which 
contribute to two- loop order. Note that there is a sum- 
mation over the lead index in Eq. (I90[) . Therefore an 
electron-hole pair created in e.g. the right lead due to a 
finite voltage enters the scaling equation for the coupling 
to the left leads. 

The signatures of the transconductance can thus be 
observed as a cutoff of the divergent coupling to the right 
lead even when the voltage is only applied to the left 
lead, see Fig. [SJ Note that the decoherence scale though 
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FIG. 5: Flow equation solution using the mean value of the 
couplings Qmean and — 2p^,"' r ? to the right leads; solid lines 
illustrate the flow for voltage Vl = Vr — 0. In this case the 
singlet-triplet excitation energy K provides a cutoff scale in 
j™';J. For a finite voltage applied to the left lead, i.e. Tk = 
0.005, N(0)K = 0.025, Vr = and N(0)(eV L ) = V = 0.05, 
we show again g^ean and — 2p^™',f (dashed lines). A decay 
of the couplings sets in which is initiated by the decoherence 
on the left side of the DQD. See discussion of the transcon- 
ductance effect in the text. For comparison we also show that 
the decoherence rate for j™™'^ and 2p^™' r f' (dotted lines) is 
larger. 

is larger than the decoherence on the right side. To study 
the transition from a strong-coupling to a weak-coupling 
problem, the calculation of either the transconductance 
or another physical quantity is necessary. This is not the 
focus of this paper. Here we concentrate on a comparison 
between the flow equation method and pRG scaling. We 
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have shown that they are identical to lowest order and by 
Fig. [5] also that the flow equation method provides the 
same physics as was found from previous studies of the 
double quantum dot system using pRG 14 . 



C. Comparison 

As a further comparison we show the mean values of 
the couplings in the two different calculations. The aver- 
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Comparison of the averaged coupling functions in 
tt]mean and g s t =2 [g s t ] mean ) and 
9mI2n and p=P™ ea „), same param- 



FIG. 6 

pRG (solid lines: g u 
flow (dashed lines: g sum 
eters as in Fig. [3] and [4] Inset: comparison of the different 
cutoff functions, X and exp(— (x/A) 2 ) 

aged quantities enter physical quantities like the current 
and spin susceptibilities. They are defined as 



sumj 
hfmean 



1 



Vj + 2K J_ 



Vj/2+K 



mean t/ ( I 

Vj + A J-Vj/2-K/S 



Vj/2-K 
V./2+K/2 



j su 

de k g k 



dc 



k Pi 



(96) 
(97) 



and analogous for the mean values of g^ t m (uj) and <7™ t m (u;). 

In Fig. [6] we show a comparison of the triplet-triplet 
and singlet-triplet couplings versus the cutoff /flow pa- 
rameter for the same set of parameters as discussed be- 
fore. We observe that the couplings in both methods 
start to grow logarithmically. The quantitative behavior 
is slightly different due to the different cutoff schemes. 
Note that the plotted variable is the averaged value over 
a frequency / energy regime which thus includes the slopes 
of the most divergent coherent couplings, compare Fig. [3] 
and |U These are different due to either a cutoff of 
e -0/A) or q(a — x) as illustrated in the inset of Fig. [5] 
Both methods are limited to weak-coupling and thus the 
value of g^cm changes only by roughly 25 % and stays 
well below the strong-coupling limit. 

Note that we compare the two-loop flow equations with 
the decoherence-cutoff corrected pRG. The philosophies 
of the two different methods are obvious in Fig. [5] In 



the pRG the flow continues to grow logarithmically and 
as soon as the reduced band reaches the cutoff, the flow 
stops and the coupling stays constant. This is the value 
which is then inserted into the expression of physical ob- 
scrvables like the current or the transconductance, etc. 
In the flow equation method the decoherence enters dif- 
ferently. As soon as the flow parameter B reaches 1/T 2 , 
the initial couplings start to decay again. Once the de- 
coherence scale is reached the Hamiltonian changes its 
form, i.e. in the Kondo model the impurity spin is en- 
tangled with the leads and therefore the newly generated 

(2) 

couplings in fl> n j start to grow and determine the dy- 
namics of the system. For B — > oo the Kondo coupling 
would thus flow to and the physical relevant value of 
the Kondo coupling should be chosen as the maximum 
before the value starts to decay. 

The two decoherence scales of the two different meth- 
ods for the example in Fig. [6] seem very different since 
the pRG flow still continues while the couplings in the 
flow equation methods already start to decay. This ob- 
servation can be traced back to two reasons. First, it is 
not obvious if B = 1 /A 2 is identically fulfilled or if there 
is some prefactor involved which invalidates the direct 
comparison. Second, as pointed out before an additional 
cutoff is found in the flow equation method which is pro- 
portional to the exchange interaction K and not given 
by the noise fluctuations. This term is not found from 
a Korringa-rate calculation in the pRG method. On the 
other hand the term contains the physics of the quan- 
tum phase transition in the double quantum dot system 
since the Kondo coupling is expected to diverge only if a 
Kondo singlet is built up with the conduction electrons 
in the leads, and not if the quantum dot is locked in a 
singlet configuration. Detailed study of the behavior of 
this transition is the subject of a future publication^!. 

The best comparison which is thus possible between 
the two methods is the frequency-dependent coupling 
directly at the decoherence scale in the flow equation 
method and for A — > in the pRG approach. This is 
shown in Fig. [7] The parameters are chosen such that 
the voltage is the determining cutoff scale and are iden- 
tical to the parameters in all other figures. Note that we 
only show the coupling in the lead where the voltage is 
applied in contrast to the discussion of the transconduc- 
tance. We can conclude that the flow equation method 
and the pRG approach describe the same physics of de- 
coherence contributions out of equilibrium. 



VII. CONCLUSION AND OUTLOOK 

In conclusion we have shown for the example of the 
double quantum dot system that the flow equation 
method^ and the poor man's scaling approach to non- 
equilibrium 2 " are equivalent to one-loop order. In both 
methods we find that the coupling develops a frequency 
dependence and that only at energies / frequencies where 
coherent processes are possible a logarithmic divergence 
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FIG. 7: Comparison of the coupling g\ t (u>) (gl}(u)) in the 
pRG (solid lines) for A — > and g k UTa ' (Pk) for the flow 
equation (dashed lines) at the decoherence scale B — 1/F 2 . 
Parameters are chosen to be the same as in the previous fig- 
ures. Note that we compare a 2-loop calculation (flow equa- 
tion) to an effective theory which mimics the effects of 2-loop 
contributions. 

of the coupling can be observed. 

If the infinite series of coherent processes is broken by 
decoherence, e.g. like in the presence of a not negligible 
current, the divergence should be lifted. The two dif- 
ferent methods use different approaches for this: while 
the pRG approach includes a physically motivated cutoff 
given by lifetime broadening and/or vertex corrections, 
the flow equation method is continued to higher order, 
i.e. two-loop, where a cutoff arises systematically. 

The double quantum dot setup is interesting to study 
since it is the simplest model which includes spin- 
coherent inelastic scattering processes like the singlet- 
to-triplet transition at the exchange energy K. The 
triplet-triplet coupling diverges at the Fermi energy if the 
triplets are degenerate, whereas the singlet is the non- 
degenerate ground state if K > 0. The lifetime of the 
ground state is infinite if there is no current applied and 
thus the double quantum dot system is also a good case 



study for the effect of decoherence due to a finite cur- 
rent. If the current and therefore decoherence exceed a 
threshold as discussed in detail in Ref. [H, this leads to 
a non-equilibrium induced current through the exchange 
coupled quantum dot even in the linear response regime. 

A further aspect of the model is the quantum phase 
transition inherent in a double quantum system due to 
the competition between the exchange coupling induced 
singlet ground state and the dynamic singlet state due 
to the Kondo interaction with the leads. For two impu- 
rities embedded in a metal a quantum phase transition 
occurs for K w 2.2T^-, which is in a regime where the 
poor man's scaling approach breaks down. This quan- 
tum phase transition can be studied in more detail in 
the flow equation approach since i) the two-loop equation 
([Mf contains the competition between singlet formation 
of the two dots vs. individual Kondo screening and ii) the 
Kondo coupling decreases below the decoherence scale in 
non-equilibrium. In order to do so a physical quantity 
like the spin susceptibility will be studied in a different 
publication 27 . Such a calculation is based on the key fea- 
ture of the flow equation approach that decoherence and 
interaction effects, e.g. the spin-spin interaction in the 
double-dot system, are dealt with on the same footing. 

We would like to thank J. Paaske, P. Wolfle and 
H. Schoeller for useful discussions and input at various 
stages of this project. Furthermore we would like to 
thank Ch. Bruder, B. Braunecker, K. Flensberg, S. An- 
dergassen, D. Schuricht and L. Fritz for valuable discus- 
sions. 

S. K. acknowledges support through SFB 631 of the 
Deutsche Forschungsgemeinschaft (DFG), the Center for 
Nanoscience (CeNS) Munich, and the German Excellence 
Initiative via the Nanosystems Initiative Munich (NIM). 



Appendix A: One-loop expressions 

With the expression for r\ in Eq. (|3"T1) and inserting it 
into the flow equation (|35[) . we find the scaling equation 
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A T sum '3(B) 



dB 



dB 



- \ £((i - »(«')) - »(«')) ((<* ^) ^ ^)) -nT'-'T-- 1 

2(1 - rc(gj)) (fafc' - + *0 - (e, - £fc - *0) 2i^,,mj fc 
+ 2 -* q -K)~ - ^ + *)) 2m£,,i* (Al) 



-{e k ,-e k + KfPl, 



k'k 



~ 2(1 - "(«?)) ((<*< - e,) - (e, - e fe + if)) ^'^ fc 

+ 2 "(«?) - 6, + tf) - (e, - e fe )) J~ (A2) 



With the diagonal parametrization as defined in the We assume the exponential dependence on the en- 
main text we can write ergy is stronger than the energy dependence of the cou- 
pling function. Doing the integration assuming that 

dJ S k Um J (B) = e -2B(x-c ff^ „ e -2B(x~cf y( c ) and neglecting contri _ 

butions from Bq = e~ SA ° we find 

- 2(1 - 2nto))(e fe - e g )e-^-^ 2 (J^ /2 ) 2 



- 42(1 - n(qj))(e k - e q + K)e~^-^ (P( k+q)/2 f 
q 

+ 42-to)(^-^-^)c- 2B( -- e -^ 2 (pr; +9)/2 ) 2 

(A3) 

dPj{B) 
dB 

- £(! - n{qj)) (2(e fc - e, - if/2)) c -2B( efc - e ,-K/2) 2 

,sumj pj 
X J (k-K/2+q)/2^(q+k+K/2)/2 

+ £ (2(c fc - e 9 + if/2)) e^ B ( e -^ + ^ 2 ) 2 

X r (k-K/2+q)/2 J (q+k+K /2)/2 \ Aq: ) 
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dgT^jB) _ J_ / rj 2_B( £fc -\/,/2) 2 , 1 -2B(e fc + \/,/2) 2 A / sum,jN2 

, O J- ( r 3 ~2B(e k -V 1 /2+K) 2 , r 3 -2B(e fc -V, /2+if) 2 \ / j \2 

I o 1 ( rj c -2B(e t -Vj2-K) 2 , 1 -2B(e fc + ^ /2-g) 2 \ / j \2 ^ 

d Pj( B ) _ 1 ( r j -2B(e k -V,/2-K/2) 2 , 1 -2B(e t /2-K/2) 2 \ sum, j i 

dB ~2B [ l + rj + l + rj )9 k -K/7Pk 

, 1 / r J -2B( efc -y i /2+K/2) 2 , 1 -2B( efc + y,/2+g/2) 2 \ j sumj 

+ 2B [ l + rj + l + rj )Pk9 k+K /2 



(A6) 



where we have now introduced the dimensionless cou- 
plings = N(0)J s »r and p{, k = N(0)P^. This 

expression is for rj = 1 discussed in detail in the main 
text. 



be neglected. The integration over momentum yields 

J de k , J de k [n(kj)(l-n(k'j))] 

x (ey -e k + aK)c- 2B ^'- tk+aK ^ 
(l - erf{V2BaK)^ 



2ir 
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Appendix B: Flow of the RKKY interaction 

In second order in the Kondo coupling we also generate 
an RKKYdikc spin exchange interaction, Irkky SlSr. 
Together with the flow of Irkky 

-^Irkky =£ {Hk'j) n(fcj)]^ J™ 



k'k 



- \n(k'j) - n(kj)} 

x(2 ^M k ,+2^iPi k ,) 
+ [n(kj)(l - n(k'j)) + n(k'j)(l - n(kjj)} 

x ^^ k Mi k ,-A^iPi k ,)} (Bl) 



there is also a constant proportional to 3/4?i 2 generated 
in the flow to second order. 



J de k > J de k [n(k'j)(l - n{kj))] 

x (ey - e k + aK)c- 2B ^'~ €k+aK ^ 

^-l^( 1 + Crf( ^^) 
Utilizing these results we thus find 



_d_ 

dB 



Irkky [l - ^(V2BK)} {2f max f (B3) 

-GC£ J ') 2 } (B4) 



dB 4 



fi 2 ^{[n(fc',)--(fcj)]<T J ^ m 



k'k 



d tp _ 3fc2 V27T f sum,n2 (r, j \2\ 

dB ° nSt ~ 4 8B 3 / 2 I ^max J K^Pmax) j 

(B5) 

where g^^x 1 an d Pmax are the couplings at the most di- 
vergent energy argument. 

Note that for BK 2 > 1 where eri(V2BK) w 1, the 
following relation holds 



+ [n(k'j) n(kj)](2 rfJ k M{ kl + 2 i^J*,)} 

(B2) 



In the following we study only the case for zero voltages, 
V L = Vr = 0. 

If we now insert the diagonal parametrization for the 
couplings and do the energy integration assuming like 
before that the energy dependence of the couplings can 



3 . 2 dlRKKY dE const _ 

4 dB dB ~ 



(B6) 



i.e. the combination of the two newly generated couplings 
stops to flow as soon as B has reached the energy scale 
l/K 2 . Corrections to each of the couplings Irkky and 
Econst are given by a term proportional to B~ z l 2 which is 
negligible small compared to the logarithmic divergence 
for small values of B. 
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Furthermore in the opposite limit BK 2 <C 1 



dl 



RKKY 



dB 



^{ [ - 1 + ^V2BK] (2p 3 max ) 2 + (g^ J ) 2 } 



K 

- — (On 3 \ 

D \ rmaxl 



(B7) 



where we neglect the small term cx 1/B 3 / 2 . 

Note that the frequency dependence of the coupling 
functions is not developed in the limit BK 2 ^ 1 and 
thus we find that the relation g^™J = 2p 3 max is fulfilled 
during the flow. Thus for example the coupling g^^ x 3 
fulfills the flow equation 



where the factor of 2 originates from the summation over 
the two lead indices j which is special to the chosen 
model. The latter equation can be solved immediately 
and yields 



Irkky(B) = K 1 



H^BTk) H^B^Tk] 



(Bll) 



where the initial value is given by Irkky {Bo) = K and 
Bq = I/Aq. The above sketched calculation is only valid 
for BK 2 <C 1 but as we have also argued the flow for 
BK 2 ^> 1 is negligible and thus we can write 



f ] a sum,j 
zniax 

dB 



— / sumJ ) 2 



sumj _ 
Umax 



HBT 2 ) 



(B8) 
(B9) 



where Tk is the Kondo temperature as defined in the 
main text. In the following we assume that the initial 
value K on the right hand side of the flow equation is 
also subject of the flow and thus replace: K — > Irkky- 
Consequently we have to solve the differential equation 



d In Irkky ~ -2 



1 



hbt 2 : 



dlnB 



(BIO) 



Irkky(B -> oo) « K[l - { 



MK/Tk) 



ln(A /T K ) 

(B12) 



The calculation in this appendix illustrates that there 
are terms which lead to a logarithmic correction of the 
exchange energy gap. The renormalization of the quan- 
tum dot energy levels is known and occurs for example as 
the Knight shift for a spin-1/2 quantum dot in magnetic 
field!. 

Note that there is no further contribution to Irkky 
from the two-loop calculation. 
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Appendix C: two-loop 



(2) 

The new interaction terms in H \' t as defined in Eq 



are generated due to two Kondo spin scattering processes 
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dB 
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dB 



= -(e k > -e k + e q , - e q - K) (K k , k . qlq 

+ ((e k ,-e k )-(e q ,-e q -K))Jl™'iM g !, q 
+ ((e k ,-e k -K)-(e q ,-e q ))Mi, k J^ 



K M N 
^q'q; k ' k ) 



(C3) 
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J( K RKKY , K RKKY\ 

"■y^k'kw'q T q' q;k' k J _ , \2 / t^RKKY ■ tsRKKY\ 
— —\ e k' - £fe "I" V — EgJ {J^k'k;q'q + ^ q' q-k' k ) 

- 2((e fe , - e fc + A") - ( V - 6, - K))P^M v q , q 

+ 2((e k , -e k -K)- (e q , - e q + K))M 3 klk P^ q (C4) 
Exchange of summation indices provides us with a symmetry constraint 

t-snm _ I/sum (OK} 
^k'k-.q'q — ^q'q-.k'k' \^°) 

rP/M _ _j^P/M 
v k'k\q'q q' q\k' k^ 

RKKY _ T^RKKY 
^k'k;q'q ^q'q;k'k 

In addition the hermiticity of the Hamiltonian has to be fulfilled and such there are some simplifying relations 

^k'k:,q'q — -^kk';qq' {^°) 
^P/M _ jy-M/P 
v k'k;q'q kk' ;qq 

RKKY _ tsRKK 
v k'k;q'q ~ ^kk'-qq 



K Jk;q'q ~ ~ K q'{;k'k' (^6) 

t^RKKY _ RKKY fr ,~\ 
^k'k-q'q — ^q'q-.k'k ■ ' ) 



K P/M _ _ R M/P (CQ , 
^k'k-q'q — ^kk'-qq' I 

K k'k3q = K k^Mq^ (CIO) 



These relations are also fulfilled in the scaling equations. 

The canonical generator r/ 2 ) = [Ho, Hf nt ] in second order of the Kondo coupling is explicitly given by 

Vtnt =( e k> - £k + Eg' - Eg)^fc™g'g* '■ (^L + Sr) ($(k'j)(kj) X «(<?'«)<») ■ 

+ (Cfc' - £fc + Eg' - £g + K)Ky k . q , q i \ (^(S L - Sr) + 2i(S L X Sr)^J (s(fc'j)(/y) X S(q'v)(qv)) : 

+ (efe' - £fc + e g ' - £g - K)Ky k . qlq i : {{^Sl - Sr) - 2i(§L x Sr)^ (&(k'j)(kj) x s*(g'„)(g„)) : 

+ (efe' - £fe + £g' - e q)K$ k .p q i2i : (S L x S(k'j)(kj)) (Sr x + (<Sfl x S(k'j)(kj)){SL x S(g'„)(g„)) : 

(Cll) 

= ^7 + < 2) + + < 2 f ^ (C12) 

Note that both X sum and K RKKY do not involve a singlet-triplet transition like K p l M . 

From [r/ 2 ),ijW] and [r/ 1 ), i^ 2 '] we find the higher order contributions to the flow dH/dB and thus a scaling 
equation for the coupling to two-loop order 



sum, j 



dB 

+ \Y, W^X 1 ™M) + n(qv)(l n(q'v))} {ey e k + e q , - e q - (e q e q ,)) (k^ qlq - k^ k , k )g q ^' v 

q'q 

+ \ £(n(<z'«) - n{qv)) (e k , e k + e q , - e g - (e q e q ,)) {k^, Y q + k™^)g~ 

q'q 

+ 2 £ n(qv)(l n{q'v)) {e k , e k + e q , e q K (e q e q , + K)) {k^ q , q - k^ k , k )p v qq , 

q'q 

+ 2^2n(q'v){l-n(qv))(e k , - e k + e q > - e q + K - (e q - e q > - K)) {k% k . q , q - k p , q . k , k )m v qq , (C13) 

q'q 

where k k i k - q i q — N(0)K k > k;q > q and we neglected a contribution from 

-2 J2(n(q'v) - n(qv)) (e k , - e k + e q , - e q - (e q - e q ,)) (k^J + k^)g s q ^ : (S L + S R )s ik , j){kj) : (C14) 

q'q 

It can be shown in the simplest limit that this term is proportional to B~ 3 / 2 and thus negligible compared to the 
leading order B^ 1 ^ 2 . In the more general cases it can be shown numerically that the contribution from k k > q . q ' k does 
not have an effect on the flow. 
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Furthermore the expression for p 3 k , k yields 
dB 

+E n W( 1 

q'q 
q'q 
q'q 

where we neglect contributions from 

- 5>(</«)(l - n(qv)) + n(qv)(l n(q'v))] ((r,™ + ^D^V - + K^MX) ( C16 ) 

q'q 

x : ((5 L - S R ) + 2i(S L x S R ))s^ k ,j^ kj) : 
+ 2 ^ ^)(l - #)) ((« + V§£§m ~ {Kfiff + Kftjfi)r&) (C17) 

q'q 

x : ((& - S R ) + 2i(S L x S R ))s [k , j){kj) : (C18) 

This is a well-controlled approximation using the argument before that the terms proportional to k'q; q'k can be shown 
numerically to be negligible small compared to the terms of type k'k;q'q. Additionally we can use the symmetry of 
the system that p k , k = —p§i k and thus the first term in the latter equation cancels and can therefore be neglected. 

Since we are not interested in the flow of the newly generated couplings, we integrate out the scaling equations for 
them in order to find their B-dependence. We illustrate this procedure on Kfjj^, as an example: 

B 

(KTq'q K^ k , k ){B) = e -*(^-«*+V-O a J dB > e +B> ( e kl -e k+ e q ,-e q r 

{((e k > - e k ) ( V - e q ))j~{B')^{B') 

+ 2((e* -e k + K)- (e q , - e q - K)) P^ k (B')M v qlq (B') 

+ 2((e* -e k -K)- (e q , - e q + K))Ml, k (B')P;, q (B')} (C19) 

In the following we apply the diagonal parametrization for the Kondo couplings and thus find 

B 

(KlT qlq ~ K^ k , k ){B) = C -^ k '~e k+ e q ,-e q f J ^ &+B ' (e k , -e k+ e q , -e t )' 

- e fc ) - (e q , - e q ))e' B '^'-^ 2 Jl um '\B')e- B '^'-^ 2 J™(B') 

+ 2({e k , -6 k + K)- (e,/ - C q - K )) e -B'(e k ,-e k+ Kf p j {Bl)e -B'(e q ,-e q -Kf pv ^ 

+ 2({e k , -e k -K)- ( V - e q + K))^ 3 ' ^'-^ K ^ ^'^'(v-^) 2 P*(B')} 

(C20) 

In the calculation to one-loop order we found that the Kondo couplings depend logarithmically on the flow parameter 
and thus have a much slower dependence than the exponential function in the latter expression. In the following we 
therefore assume that we can replace the coupling by their averaged value in the interval from Bo to B. In general the 
initial value B = 1/Aq should be chosen as the band cutoff in order to treat boundary terms. Since this is equivalent 
to treating a different model with new boundary conditions we set in the following as also in the main text 

B = (C21) 



- n(q'v)) (e q , - e q + e k> - e k + K - (e q - e q ,)) (k£ k . qlq - k^, q . k , k )g s qq f 1 ' v 

- n(qv)) {e q , -e q + e k , - e k - (e, - e q , + K)) (k s k ^ qlq - k s q Y™ k , k )p v qq , 

- n{qv)) (e q , - e q + e k , - e k - (e, - e q , + K)) (fc™ Y + k^l)p v qq , (C15) 
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Using these two simplifications the integration is simple and can be done straightforwardly yielding 



i.K^ qlq -K^ k , k ){B) = _L e - B < £ *'- £ *+V-«.) 



(efe' - £fe) - (e q ' - e q ) 



{ l e fc' - e 
I (ey- 



«L_!iZ(e 2B (^'-«)(v-^)-i) ± dB'J s k um ^(B')J™ m > v (B') 
(ey - efe)(e 9 ' — e 9 j LB 



, (efc' - e fc + AT) - (e g / — e q — K) 
(efe' - £fc + ^ )(V ~e q - K) 



( e 2B(t k ,-e k + K)(e q ,-e q -K) _ ^ 



- / dB'Pl{B')pv{B') 



+ 2 (efc' - e fc - AT) - (y - £ g + jQ , 2B (efc ,- e . -jQ( v - e .+«:) 
(e fc , - e 9 + X) ^ 



JdB'P>(B')PZ(B')]} 
(C22) 



The same procedure is applied to the other coupling generated to two-loop order and these expressions are then 
inserted in the flow equation for the Kondo coupling. In the following we write down only the result for the diagonal 
parametrization where we assume e.g. in g s k umj that ey = € k . We find that the collected terms contributing to 
dg s k um ' 3 jdB in two-loop order can be written by TL(q, q') + TL(q\ q) where 



TL(q,q') = ^^n(^)(l - n(qv)){ - 2(e q , e q ) 2 ) C - 2B ^"^ 2 

q'q 

[ f 7 r~il sum. 7 / t-)/\ sum.f / ~r>! \ sum.f / t-)\ 
- dBg k "(B)g (ql ; q)/2 {B)\g {q+ql)/2 {B) 

+ 2Bj2n(q'v)(l - n(qv))( - 2( V - e q + K f) c -^(, q , -e q+K f 

q'q 

[1 f ^T j (b')4 + „ /2 (b')]pU)/ 2 ( b ) 

J Bq 

+ 2 £ [n(q'v)(l - n(qv))]e-^-^ ^=gj^ ~f (e 2B 



(-K)(e ql -c q +K) 



1) 



[1 f ' dBy k (B>) P J ql+q)/2 (B>) 



sum,D / p\ 



+ 2]T [n(q'v)(l - n(qv))]c- 2B ^'-^ 2 K2 e f ^B( + K) { e q ,-, q ) _ 1} 

q'q H H 

i / B ^v fe (5')v+;)/ 2 ( 5 ')]p( 9+90 / 2 (5) 

' ' ' Bq 



(C23) 



Note that TL(q, q') = TL(q\ q). For the symmetric setup we are studying, only the first two terms are non-zero after 
summation over the lead index (p k = —p k ) and the last two terms are thus neglected in the following. Otherwise the 
last two terms are nasty whereas the first two terms can be summed over momentum straightforwardly. 

Please note that the last two terms are of the same order than the first two terms if K = 0. In the expressions 
which are discussed in the main text it is thus not valid to set K = 0, especially in the case when one wants to derive 
the spin- 1/2 limiting case where the cancellation of p k = —p k does not take place. 



In dp/dB the terms are less symmetric 
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dB 



= one-loop 



+ BJ2 n(qv)(l - n(q'v)) (-2(e q , - e q f) e~ 2B ^'~^ 2 

q'q 

I ^ f Ji~>l i ( t»I\ sum.!? / t->/\1 sum.!? 

Lb L dB ^( B )9 {q '+ q )/2( B )\9 {q +^/2 

+ 4B^n(g'v)(l - n(qv))( - 2(e q , - e q - K )^ c -2B(e q ,-e q -Kf 



q'q 



P(q+q')/2 



[1 f dB>p> k (B')pl ql+q)/2 (B>) 
+ ^n(q'v)(l - n(^))c- 2B (v-^-^) 2 l!«_-iill^( c 2B(-^)(v-^) _ i) 



B 
Bo 



+ ^2n(gv){l - n{q'v))c- 2B( ^'- e ^ 



P{q+q')/2 



q'q 



9 s r^B') P \ ql+q)/2 {B') 



(e q , -e q + K) 2 - K 2 
(e q '- 

sum.D 



K(e q , -e q + K) 



^2B(-K)(e g ,-e q +K) _ ^ 



'(9+9')/2 



(C24) 



Note that we do not have q o q' symmetry! This is 
somehow expected since pk'k is only symmetric under 
exchange of k <H> k' and also K o —K. The latter 
symmetry is fulfilled in the scaling equation. 

For the energy integration the following relation is use- 
ful 

J de q , J de q n(qv)(l - n(q'v)) (C25) 

x (e q , -e q + aK) 2 c- 2B (^-^ +aK ) 2 
= il-^_)F {a K) (C26) 



+ 



(i + r v y 



(F(aK + eV v ) + F{aK - eV v )) 



(C27) 



where contributions from the band cutoff cancel out or 
are exponentially small and 



f{v) = 



y 



8BJ2B 



erf 



(V2By) 



+ 



2(2B) 2 



-2By 2 



(C28) 



The explicit expression for dg^™ 1 ' 3 /dB after integration 
is given in the main text. 
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and 



where 



7 sum, j / ni 

a 9k \ U ) _ 1 ( r 3 -2B(e k -Vi/2) 2 , 1 -2B(e k +V, /2) 2 \ / sum,j\ 



2 



dpj(B) 
dB 



2BW + rj 

2 1 ( T 3 c -2B(t k -V,/2+K) 2 + r J c -2B(£ fc -y,/2+A')- 



2B V 1 + 1 + rj 



o 1 f ^ r -2B( £fc -V,/2-K) 2 1 -2B(c k + V 1 /2-K) 2 \ ( j \" 



-E ^r ,j ( 5 )te m /(s)) 2 

V 

9(0,r v ,V v )g s r' j (B)( 9 s ™> v (B)) 2 

V 

E (^~ 2BK2 + f T^erf (^/f))^(B)(2pS.« (S)) 2 
- £ ff(if,r„,K) 5 r iJ (S)(2^a,(^)) 2 (C29) 

V 

- 1 f r 3 -2B(e k -V 1 /2-K/2) 2 , 1 -2B( efc +V, /2-K/2) 2 N sum.j j 

■2B l l + ri + l + rj J^-*/2^ 

I 1 / r J -2B( et -y i /2+K/2) 2 , 1 ~2B(e k + V n /2+K/2) 2 \ j sum.j 

+ 2i?U + r J + l+r,- ^fe-W/2 

-E ^i-( B )(5~(s)) 2 

- ^ g(0, r„, V v )p{(B)(g°™f(B)) 2 

V 
V 

-J2 g(K,r Vl V v )p{(B)(2 P : nax (B)) 2 (C30) 



g(K,r,V) = ^^_(e- 2S (^) 2 + c « + y) 2 _ 2e -2^ 2 } 
4i> (1 + r) J v ' 

' ' " ' :{{K- V)etf(V2B(K - V)) + (K + V)erf (Vw(K + V)) - 2Ker{(V2BK)) (C31) 
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For r = 1 the function g(if , V") corresponds to the function f(K, V) in the main text. 
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